A = sym([1, 3; 2, 2])
[P, D] = symeig(A)